Head-on collisions of unequal mass black holes in D = 5 dimensions 
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We study head-on collisions of unequal mass black hole binaries in D = 5 spacetime dimensions, 
with mass ratios between 1:1 and 1:4. Information about gravitational radiation is extracted by 
using the Kodama-Ishibashi gauge-invariant formalism and details of the apparent horizon of the 
final black hole. We present waveforms, total integrated energy and momentum for this process. 
Our results show surprisingly good agreement, within 5 % or less, with those extrapolated from 
linearized, point-particle calculations. Our results also show that consistency with the area theorem 
bound requires that the same process in a large number of spacetime dimensions must display new 
features. 
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I. INTRODUCTION 

Black holes (BHs) have been at the center stage of fun- 
damental physics in the last decades: supcrmassive, as- 
trophysical BHs lurk at the center of most galaxies while 
large numbers of stellar-mass BHs are thought to popu- 
late each galaxy [l|, Q ; highly dynamical BH binaries are 
strong sources of gravitational waves and, perhaps, power 
jets and other extreme phenomena 0, Q • I n high-energy 
physics, BHs are a central piece of the gauge-gravity du- 
ality @, and arc the generic outcome of particle collisions 
at center-of-mass energies above the Planck scale @ . In 
this regime the particular nature of the particles' struc- 
ture should become irrelevant, as indicated by Thome's 
"hoop" conjecture Q, and "no- hair theorem" -type of ar- 
guments. These can be invoked to argue that, in general, 
trans-Planckian collisions of particles are well described 
by collisions of highly boosted BHs. In this context, sce- 
narios such as TeV-gravity are especially interesting, as 
they lower the Planck scale to the level at which BHs 
would beproduced in cosmic rays and particle acceler- 
ators |8l-[r7|. Thus, high-energy BH collisions could be 
used to look for signatures of extra dimensions and BH 
production in ground-based experiments in the forthcom- 
ing years. At the fundamental level, BHs might hold the 
key for a theory of quantum gravity, and might help un- 
derstand important issues such as cosmic censorship, in- 
formation loss and the maximum possible luminosities in 
any physical process [HI, [l9[ . 
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The above arguments illustrate the necessity to un- 
derstand accurately dynamical BH spacetimes, and their 
potential across a wide variety of fields. Since the full 
system of Einstein equations needs to be carefully under- 
stood, this is by all means a monumental task, and typi- 
cally requires numerical methods. With these fundamen- 
tal issues as motivation, long-term efforts to understand 
dynamical BHs in generic spacetimes have been initi- 
ate d l20l - [2o| , ranging from the inspiralling of BH binaries 
ll-Sir high-energy collisions of BHs in four 0,[ll[2| 
and low energy collisions in higher spacetime dimensions 
[III, [HI, stability studies in higher dimensions [30l - [32| and 
BH evolutions in non asymptotically flat spacetimes [33j . 

Our group has recently studied head-on collisions of 
equal-mass black holes in higher dimensions, in partic- 
ular D = 5 [U [22| , hereafter denoted as Paper I and 
Paper II, respectively. In the present work, we wish to 
extend that study to the case of unequal mass BH bina- 
ries. This is an interesting extension for several reasons, 
perhaps the most important of which is the nontrivial 
comparison with point-particle (PP) calculations in the 
linearized regime. We will compare radiated energy, mo- 
mentum and multipolar dependence of our full nonlinear 
results with results from linearized Einstein equations. It 
turns out that the agreement is remarkable, providing an 
outstanding consistency check on our codes and results. 
A thorough analysis of the linearized Einstein equations 
is done in an accompanying paper (34j . 

This paper is organized as follows: In Sec. |TT]we sum- 
marize our numerical method and setup and present the 
numerical results. We finish by giving some conclusions 
and final remarks in Sec. Mil 
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II. NUMERICAL RESULTS 



The numerical simulations have been performed with 
the Lean code, originally introduced in Refs. [H, HU 
and adapted to higher dimensional spacctimcs in Pa- 
per I. The Lean code is based on the Cactus computa- 
tional toolkit |37l| and uses the Carpet mesh refinement 
package [38, 39] and Thornburg's apparent horizon finder 
AHFinderDirect jlfllZll]. 

Following the approach developed in Papers I and II, 
we perform a dimensional reduction by isometry of a D- 
dimensional {D > 5), asymptotically fiat spacctime with 
SO(D — 2) isometry group; this is the symmetry, for in- 
stance, of a head-on collision of black holes in D dimen- 
sions. As discussed in Paper I, the dimensional reduc- 
tion is performed on the (D — 4)-sphcre S D ~ A , which 
is the symmetry manifold generated by the subgroup 
SO{D - 3) C SO(D - 2). We remark that the dimen- 
sional reduction is not a compactification, but simply a 
way to employ the symmetries of the problem in order to 
rewrite the D dimensional vacuum Einstein equations as 
an effective 3+1 dimensional time evolution problem with 
source terms that involve a scalar field (see Eqs. (2.16)- 
(2.18) in Paper I). 

The coordinate frame in which the numerical simula- 
tions are performed is 



^ D - i ) = {t^y,z,cp\...,4> D - i ) 



(1) 



where the angles cj) 1 , . . . , <f> D ~ A describe the quotient man- 
ifold S D ~ A and do not appear explicitly in the simula- 
tions. Here, z is the symmetry axis, i.e. the collision line. 

We have evolved this system using the Baumgarte- 
Shapiro-Shibata-Nakamura [42|, |43| formulation along the 
lines presented in Paper I (see eqs. (2.43a)-(2.49b) 
therein), together with the moving puncture approach 
flil |45| . The initial data consist in the time-symmetric 
Brill-Lindquist initial data in the form presented in Pa- 
per II (Eq. (2.15) therein). Gravitational waves have 
been extracted using the Kodama-Ishibashi (KI) formal- 
ism (4(| |4?| . For details of the wave extraction implemen- 
tation we refer the reader to Paper II. We have evolved 
BH binaries, colliding head-on from rest with mass ra- 
tios q = Mi/M 2 = rgi 3 / r s,2 3 = !> V 2 > V3, 1/4, where 
Mi is the mass of the i-th BH. The mass parameter 
r S,i 3 / r s~ 3 of the smaller BH is given in Table U and 
we adapt the value of the second BH accordingly. The 
initial coordinate separation of the two BHs is set to 
d/rs = 6.37 which translates to a proper initial sepa- 
ration of L/rs = 6.33. Further details of the setup of the 
simulations are summarized in Table |TJ Unless denoted 
otherwise, our discussion will always refer to the highest 
resolution runs with hf/rs = 1/84, hf/r$ = 1/102.9, 
h f /r s = 1/118.8 and h//r s = 1/132.8 for models HD5a, 
HD5b, HD5c and HD5d/ in Table H respectively. The 
energy flux is computed according to Eq. (2.56) in Pa- 
per II (see Eq. (21) in Ref. [48[ for the corresponding 
expression in Fourier space). The momentum flux can 
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FIG. 1. (Color online) From top to bottom: I — 2, I = 3 and 
I = 4 modes of the KI waveform for the different mass ratios; 
q = 1 (black solid lines), q = 1/2 (red dashed lines) and q = 
1/3 (green dash-dotted lines). The curves have been shifted in 
time such that the formation of the common apparent horizon 
corresponds to At/rs = and taking into account the time 
that it takes the waves to propagate to the extraction radius. 



be obtained from 



dP l 



jn d2E i 
dVL n 

dtdVL 



(2) 



with n % a unit radial vector on the sphere at infinity . 
This results in an infinite series coupling different mul- 
tipoles. Using only the first two terms in the series, we 
find, for instance, that in D = 5 the momentum flux in 
the collision direction is given by 



dP 
~dt 



1 

47T 



$|r 3 ( 5$ !r 2 + 2i$ !r 4 



(3) 



Here, is the /—pole component of the KI gauge- 
invariant wavefunction [22T |46| . |47| . From the momen- 
tum radiated, the recoil velocity of the system can be 
obtained as 



00 , dP 
dt— 

.™ at 



(4) 



A. Waveforms 

In Fig. Q] we show the I = 2,3,4 waveforms for dif- 
ferent mass ratios, zoomed in around the time of the 
merger. The waveforms have been shifted in time such 
that At/rs = (t — r ex — tcAH)/rs = corresponds to 
the time tcAH at which the common apparent horizon 
forms and taking into account the propagation time of the 
waves to the extraction radius r ex /rs — 60, 49, 42.4, 37.9. 
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Run 


q 


D-3 I D-3 
r S,l l r S 


D-3 

r s 


zi/rs 


Z2/rs Grid Setup 




HD5a 


1 


0.5 


1 


3.185 


-3.185 {(256, 128, 64, 32, 16, 8) x (2, 1, 0.5), h = 1/84} 


HD5b 


1/2 


0.33 


1.5 


4.247 


-2.123 {(209, 104.5,52.3,26.1, 13.1,6.5) x (1.6,0.8,0.4), h 


= 1/102.9} 


HD5c 


1/3 


0.25 


2 


4.777 


-1.592 {(181.0, 90.5,45.3,22.6, 11.3) x (2.8, 1.4,0.7,0.4), h 


= 1/118.8} 


HD5d c 


1/4 


0.2 


2.5 


5.096 


-1.274 {(161.9,80.9,40.5,20.2,10.1) x (2.5,1.3,0.6,0.3), h 


= 1/113.8} 


HD5d m 


1/4 


0.2 


2.5 


5.096 


-1.274 {(161.9,80.9,40.5,20.2,10.1) x (2.5,1.3,0.6,0.3), h 


= 1/123.3} 


HD5d/ 


1/4 


0.2 


2.5 


5.096 


-1.274 {(161.9,80.9,40.5,20.2,10.1) x (2.5,1.3,0.6,0.3), h 


= 1/132.8} 



TABLE I. Grid structure and initial parameters of the head-on collisions from rest in D — 5. The grid setup is given in terms 
of the "radii" of the individual refinement levels, as well as the resolution near the punctures h, in units of rs (see Sec. II E in 
[H| for details). We give the Schwarzschild radius r^ -3 = r§ 'j 3 + r s,2 3 of tne nna l BH , the mass of the smaller BH rgi 3 /r§~ 3 
and the quantity q = M1/M2 = fj^/rjj 3 denotes the mass ratio. Zi/rs is the initial position of i-th BH. 



The waveform is similar to previous four-dimensional re- 
sults (see, e.g. Ref. [i^] , a more detailed study is in prepa- 
ration [50] ) - Although not shown in Fig. [T] we observe a 
small, spurious signal starting around (t — r ex )/rs ~ 0, 
which is an artifact of the initial data. 

The actual physical part of the waveform is dominated 
by the merger signal at At/r$ ~ followed by the quasi- 
normal ringdown. We estimate that the different ring- 
down modes are given by 

uj l=2 r s = 0.955 ± 0.005 - i(0.255 ± 0.005) , 
aj l=3 r s = 1.60 ± 0.01 - i(0.31 ± 0.01) , 
uj l=4 r s = 2.25 ± 0.03 - i(0.35 ± 0.05) . 

These results agree well, and within uncertainties, with 
estimates from linearized theory [H, [5lTl53l |. providing 
a strong consistency check on our results. Finally, we 
consider numerical convergence of our waveforms. This 
study is summarized in Fig. [2] for the I = 2 mode of 
the KI wavefunction, and for the most challenging mass 
ratio, q = 1/4, model HD5d in Table Q] We have 
evolved this setup at three different resolutions, namely 
h c /r s = 1/113.8, h m /r s = 1/123.3 and hj/r a = 1/132.8, 
which we will refer to as "coarse" , "medium" and "high" 
resolution in the following. We show the difference be- 
tween the coarse and medium as well as between the 
medium and high resolution waveforms. The latter has 
been amplified by the factor Q = 1.47, which indicates 
fourth order convergence. We obtain the same order of 
accuracy for the higher modes. The discretization error 
in the waveforms is estimated to be sa 1.5%. 



B. Radiated energy 

Tabic [II] lists some of the most important physical 
quantities which characterize the head-on collision of BHs 
in D = 5. In particular, we show the radiated energy in 
units of total mass M, and the recoil velocity of the final 
BH in km/s. The maximum amount of energy is emit- 
ted in the equal mass case (E rad /M = 0.089% as found 
previously in Paper II) , and it decreases for smaller mass 
ratios. We estimate the error in the radiated energy to 




FIG. 2. (Color online) Convergence analysis of the I = 2 
mode for model HD5d in Table [I] We show the differences 
between the coarse and medium resolution waveform (black 
solid line) and the medium and high resolution waveform (red 
dashed line) . The latter has been amplified by the factor Q = 
1.47, indicating fourth order convergence. The curves have 
been shifted in time such that the formation of the common 
apparent horizon corresponds to At/rs — and taking into 
account the time that it takes the waves to propagate to the 
extraction radius. 



q 


E lad /M(%) 


E\l%{%) E^(%) ££ d 4 (%) 


v ieco n(km/s) 


1/1 


0.089(0.090) 


99.9 


0.0 


0.1 


0.00 


1/2 


0.073(0.067) 


97.7 


2.2 


0.1 


11.37 


1/3 


0.054(0.051) 


94.8 


4.8 


0.4 


12.64 


1/4 


0.040(0.035) 


92.4 


7.0 


0.6 


11.38 



TABLE II. Summary of our results concerning unequal mass 
head-on collisions of BHs in D — 5. We show the total radi- 
ated energy E/M as measured from the energy flux at r ex , 
the quantity in parenthesis refers to the estimate obtained 
using properties of the apparent horizon (see Paper II for de- 
tails). The next three columns show the fraction of energy 
E excited in the Z-th mode as compared to the total radiated 
energy. The last column refers to the recoil velocity 

^recoil 1H 

km/ 8. 
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FIG. 3. (Color online) Top: Total integrated energy for dif- 
ferent mass ratios, as function of if = [q/ (1 + q) 2 ] 2 . The black 
diamonds denote the numerical data and the black dashed- 
dotted line is the corresponding fitting function, Eq. ([5}. Bot- 
tom: Recoil velocity v rcco n in km /a as function of the mass 
ratio q. The red circles denote the numerical data and the 
red dashed line is the corresponding fitting function, Eq. (5). 



be about 5 %. These results have been obtained by inte- 
grating the energy flux as given by the KI master wave- 
function. We have also estimated the radiated energy 
using properties of the apparent horizon as described in 
Paper II. We estimate the discretization error to be about 
~ 10 % when using this method. The apparent horizon 
estimate for the total radiated energy is shown in paren- 
thesis in Table [ill an d is consistent with the flux com- 
putation within numerical uncertainties. Table HI] also 
shows the fraction of energy emitted in different mul- 
tipoles. Higher multipoles are clearly enhanced as the 
mass ratio decreases, in agreement with what we expect 
in the extreme case of a PP falling into a BH. In fact, we 
can make this statement more precise. Post-Newtonian 
arguments, which extend to generic D-dimensions, allow 
one to expect the functional dependence for the total ra- 
diated energy E lad /M oc if, where -q = q/(l + q) 2 is 
the dimensionless reduced mass. For clarity, we show the 
ratio E rad / (Mrj 2 ) in the top panel of Fig. [3l which can 
be seen to depend very weakly on rj 2 . We can refine this 
argument by fitting our numerical results to an improved 
expression of the form E rad /Mrj 2 = Aq + Air/ 2 . We find 



E 



rad 



Mr? 



0.0164 - 0.0336^ 



(5) 



Moreover, the following expressions for the multipolar 
content provide a good fit to our numerical data, 

Israel prad 

= 0.79 + 0.83,,, J=f =0.19-0.77*?. (6) 

Linearized, PP calculations presented in an accompany- 
ing paper [34j show that in the limit of zero mass ratio 



FIG. 4. (Color online) Momentum flux for different mass 
ratios. The curves have been shifted in time such that the 
formation of the common apparent horizon corresponds to 
At/rs = and taking into account the time that it takes the 
waves to propagate to the extraction radius. 



one obtains 



E 



rad 
PP 



Mr] 2 



= 0.0165. 



which agrees with the extrapolation of our numerical re- 
sults within less than 1%. The multipole contents in the 
PP limit are 



prad 

= 0.784, 

E pp 



rprad 

^=2 =0.167, 

tp rad 

Epp 



(8) 



still in very good agreement with the extrapolation of 
our full numerical results to the zero mass ratio limit. 
In fact, bearing in mind that we are extrapolating from 
mass ratios of 1/4 down to the zero mass ratio limit, 
the agreement is impressive. Finally, all these results are 
consistent with the fact that higher multipoles contribute 
more to the radiation than in D = 4, where for instance 
the I = 3 mode contributed roughly 10% of the total 
energy in the PP limit (5f|. Linearized, point-particle 
calculations show that the trend is consistent and con- 
tinues in higher- D [34| , which might mean that accurate 
wave extraction will become extremely difficult, as higher 
resolution is necessary to resolve higher-Z modes (56|. 



C. Radiated momentum 

For unequal-mass collisions, the asymmetric emission 
of radiation along the collision axis causes a net momen- 
tum to be carried by gravitational waves. As such, the 
final BH will "recoil", according to Eq. Q. Momen- 
tum fluxes and recoil velocity for different mass-ratios 
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are shown in Fig. [4] and the bottom panel of Fig. |3l re- 
spectively. We estimate the errors in the recoil velocity 
to be k 5 %. The general functional form for the de- 
pendence of momentum on the mass parameters of the 
individual holes has been worked out by M. Lemos in 
generic spacetime dimensions (H3 | and is the same as in 
four dimensions 



C 



g 2 (l-g) 

(i + # 



(9) 



By fitting this function to our numerical data, we obtain 
C = 716 km/s. Observe that w r ocoii reaches a maximum 
value at q = 2 — if ~ 0.38, where ip is the golden ratio. 
The quality of the fit can be seen in the bottom panel of 
Fig. [31 where we overplot the numerical data points with 
the fitting function, Eq. (|9|). This exercise is interesting 
because we can again extrapolate our results to the PP 
limit. In an accompanying paper, Berti et al. [34| find 



"recoil = 779g 2 km/s. 



(10) 



in reasonably good agreement (better than 10%) with our 
extrapolation. We note that momentum emission is given 
by a nontrivial interference between different multipoles, 
so this is a nontrivial agreement. 



III. CONCLUSIONS 

For the first time, we have evolved unequal-mass BH 
binaries in higher dimensions, by solving Einstein equa- 
tions in the full nonlinear regime. We have focused on 
head-on collisions in D = 5 spacetime dimensions. 

The gravitational waveforms are similar to the D = 4 
counterparts [HJ , and we were able to estimate the 
ringdown frequencies of the lowest multipoles. We find 
good agreement with published values for the quasi- 
normal frequencies, extracted in a linearized formalism. 
When extrapolated to the zero-mass-ratio limit, our re- 
sults agree with linearized calculations [Hf at the % level 
or better for the energy and momentum radiated, as 
well as for the multipolar dependence. This outstand- 
ing agreement is one of the main results of this work. 

Our findings, supported by linearized analysis, indi- 
cate that the higher multipoles become more important 



for larger D. This will certainly make wave extraction at 
sufficiently large D a more demanding task, since higher 
resolutions are necessary to resolve these modes. The 
momentum structure is similar to the four-dimensional 
case, it would be interesting to understand if other as- 
pects, such as the antikick, can still be interpreted in 
like-manner (57| . Finally, it would be very interesting to 
perform an exhaustive set of simulations in hig her D: our 
results, together with linearized analysis [341 ] . suggest a 
qualitative change in radiation emission for D < 12 — 13. 
In fact, this change is required by the fact that Hawking's 
area theorem forces the total amount of gravitational ra- 
diation to decrease with D, at sufficiently large D [22], [Hj]. 
Understanding the mechanism at play requires extension 
of our results to arbitrary spacetime dimensions. 
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